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It is assumed 1 3 that the radioactive decay of 44 Ti powers the infrared, optical and UV emis- 
sion of supernova remnants after the complete decay of 56 Co and 57 Co (the isotopes that 
dominated the energy balance during the first three to four years after the explosion) until the 
beginning of active interaction of the ejecta with the surrounding matter. Simulations 4 s show 
that the initial mass of 44 Ti synthesized in core-collapse supernovae is (0.02 — 2.5) x 10~ 4 
solar masses (M ). Hard X-rays and 7-rays from the decay of this 44 Ti have been unam- 
biguously observed from Cassiopeia A only™3, leading to the suggestion that the values of the 
initial mass of 44 Ti near the upper bound of the predictions occur only in exceptional cases^. 



^ ; For the remnant of supernova 1987 A^HB a n upper limit to the initial mass of 44 Ti of < 10 
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M Q has been obtained from direct X-ray observations^, and an estimate of(l — 2)xl0 
M Q has been made from infrared light curves and ultraviolet spectra by complex model- 
dependent computations^^. Here we report observations of hard X-rays from the remnant 
of supernova 1987A in the narrow band containing two direct-escape lines of 44 Ti at 67.9 and 
78.4 keV. The measured line fluxes imply that this decay provided sufficient energy to power 
the remnant at late times. We estimate that the initial mass of 44 Ti was (3.1 ± 0.8) x 10~ 4 
M , which is near the upper bound of theoretical predictions. 

The only immediate way to determine the initial mass of 44 Ti (denoted M 44 ) in a supernova 
remnant is through its emission in the direct-escape lines at energies E{ =4.1, 67.9, 78.4, 511 and 
1,157 keV. This isotope decays according to the chain 44 Ti— > 44 Sc— )- 44 Ca, with a characteristic time 
of £44 ~ 85.0 ± 0.4 years (ref. [16]). Photons in the three X-ray lines, namely 4.1, 67.9 and 78.4 
keV, are emitted with efficiencies (average numbers of photons per decay) Wi = 17.4%, 87.7% 
and 94.7% during the first stage of the decay, and photons in the 7-ray lines, namely 511 and 1,157 
keV, are produced with Wi = 188.8% and 99.9% during the second very quick (~ 5.7 h) stage. 
When the ejecta becomes transparent for hard-X-ray and 7-ray lines (typically at t > 20 years), the 
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corresponding fluxes can be estimated as 
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where d is the distance to the remnant and m p is the proton mass. Measurements^ of the 67.9-, 
78.4- and 1,157-keV lines from Cas A have yielded the value M 44 = 1.6(+0.6/- 0.3) x 10" 4 M : 
tentative detections of the 1,157-keV line from the Vela Junior remnant 17 and the 4.1-keV line 
from the Gl.9+0.3 remnanP^have been also reported, but they are more uncertain. Assuming a 
value of M 44 similar to that in Cas A, we can expect the photon fluxes in these lines from supernova 
remnant (SNR) 1987A (d ~ 50 kpc) about 23 years after the explosion to be about F 68 — 0.93i<V 8 ~ 
0.48F 68+78 ~ 0.88^1,157 ^ 0.46F 5 n ~ 3.7 x 10~ 6 photons cm" 2 s _1 . Note that this remnant was 
the first for which the X-ray and 7-ray emission from the decay of 56 Co and 57 Co was detectecP^. 

The International Gamma-Ray Astrophysics Laboratory (INTEGRAL)^ invested consider- 
able observing time in trying to detect these lines from SNR 1987 A. The first observation was 
carried out in 2003 soon after the launch of INTEGRAL. A number of important results have been 
obtained, including constraints on the luminosity of the stellar remnant, but the dedicated expo- 
sure (~ 1.5 Ms) was insufficient for the detection of the 44 Ti lines. The next set of observations 
was carried out in 2010-11 with a total exposure of ~ 4.5 Ms. Figure 1 presents the results of our 
analysis of all these data. Images of the SNR 1987A field obtained with the IBIS/ISGRI hard-X-ray 
telescope on board INTEGRAL are shown in three energy bands, 48-65, 65-82 and 82-99 keV, 
which encompass the 44 Ti lines at 67.9 and 78.4 keV. Two known sources, the black-hole binary 
LMC X- 1 and the 50-ms Crab-like pulsar PSR B0540-69, are seen in all images. Remarkably, there 
is an obvious excess at the position of SNR 1987A, but only in the energy band that includes the 44 Ti 
lines. The location of the excess coincides within 4' with that of SNR 1987A - a typical uncertainty 
for moderately bright sources (the angular resolution of IBIS is ~ 12' full- width at half-maximum 
(FWHM)). 

The signal-to-noise ratio for this excess, S/N ~ 4.1 (related to the flux measured in both 
lines at once), corresponds to a rather small probability of detecting such an excess by chance, 
P — 2 x 10~ 5 . The probability of detecting a random excess with the same S/N ratio at an 
arbitrary position in the whole 29° x 29° IBIS field of view is much higher (it can be estimated 
as mP ~ 0.43 where m ~ (29°/12') 2 ~ 2.1 x 10 4 is the number of uncorrected spots within 
the IBIS field of view; a 'spot' is the area of the sky image that can be occupied by a single 
source). In reality, the probability of detecting the 44 Ti lines by chance is even smaller than P , 
because significant excesses were simultaneously detected at the position of SNR1987A in the 
images accumulated in the narrow energy bands 62.7-73.2 and 73.2-83.7 keV. Each of these bands 
contains one of the 44 Ti lines and they are nearly independent in terms of the energy resolution of 



IBIS/ISGRI, which is 8% FWHM at ~ 70 keV. The S/N ratios measured for SNR 1987A in these 
images are respectively 3.0 and 3.1, and the corresponding chance probabilities are Pi = 1.3 x 10~ 3 
and P 2 = 9.7 x 10~ 4 . The combined probability of detecting the 44 Ti emission by chance is then 
equal to Pi x P 2 ~ 1.3 x 10~ 6 (~ 4.7a). 

The hard-X-ray spectrum obtained with IBIS/ISGRI from SNR 1987A is presented in Fig. 2. 
We fitted it using two lines of Gaussian shape with fixed centroid energies, equal widths and self- 
consistent (according to the 44 Ti emission efficiencies) normalizations. The flux measured in both 
lines at once, -F68+ 78 — (1-7 ± 0.4) x 10~ 5 photons cm~ 2 s _1 , corresponds to a value for initial 
amount of synthesized 44 Ti of M 44 ~ (3.5 ± 0.8) x 10 -4 M (the uncertainty is purely statistical). 
This estimate can be slightly decreased if we assume the presence of some underlying continuum 
in the SNR 1987 A spectrum - for example, in the form of a power law with a photon index of 2.1 (a 
Crab-like spectrum), which should describe the hard-X-ray emission well up to ~ 80 keV (ref. |23l . 
Its possible contribution at 1 keV, (8.5 ± 3.2) x 10~ 4 photons cm~ 2 s _1 keV _1 , coincides with the 
normalization of the non-thermal X-ray component in the SNR 1987 A spectrum following from the 
XMM-Newton observations in 2009 24 . With such a component (shown in Fig. 2 by a solid curve) 
the flux in the lines decreases to ^68+78 — (1-5 ± 0.4) x 10~ 5 photons cm" 2 s -1 , which may be 
translated to M 44 ~ (3.1 ± 0.8) x 10~ 4 M . 

The SPI 7-ray spectrometer on board INTEGRAL observed SNR 1987 A simultaneously with 
IBIS/ISGRI in the 511- and 1,157-keV lines of 44 Ti; however, the smaller effective area and pres- 
ence of strong background features at 511 and ~ 1, 100 keV preclude their confident detection. 
The photon spectrum measured in the 1, 106-1, 205 -keV band is presented in Fig. 3. The solid 
curve shows the 1.7a upper limit (90% confidence) for the flux in the 1,157-keV Gaussian line, 
namely Fi 157 < 2.6 x 10~ 5 photons cm~ 2 s _1 . It can be converted to the limit for the amount of 
44 Ti, M 44 < 9 x 10~ 4 M . For comparison, the dotted curve in Fig. 3 shows the Gaussian line 
corresponding to the value of M 44 measured with IBIS/ISGRI. 

The detection of emission in the 67.9- and 78.4-keV lines allows us to measure directly the 
amount of radioactive 44 Ti synthesized during the explosion, that is, M 44 ~ (3.1 ±0.8) x 10~ 4 M Q . 
This amount formally exceeds by a factor of 1.5 to 2 the estimates based on theoretical simulations 
of explosive nucleosynthesis in supernova 1987 A and on the infrared and ultraviolet data (although 
the discrepancy represents only one standard deviation). We note that 44 Ti may be produced either 
as a result of strong deviation from nuclear statistical equilibrium during the explosion (so-called 
a-rich freeze-out or as a result of incomplete Si-burning in a thin shell surrounding the region 
of main energy releaseP^. In both cases, the simulations are not straightforward and many poorly 
known parameters (for example, the mass of the collapsed core, the maximum temperature and 
density behind the shock wave, the neutron excess and the composition of the pre-supernova inte- 



rior) could affect the results. In this sense, the value of M44 that we report here may be useful for 
studying explosive nucleosynthesis and physical conditions reached during the shock- wave passage 
in the ejecta. 

We note finally that the shock interaction of the ejecta and the surrounding matter is responsi- 
ble for the synchrotron emission observed at present from SNR 1987A in the radio- and soft-X-ray 
bands 26 . The reprocessing of this emission in the ejecta is contributing to the late-time optical- 
infrared emission of SNR1987A 15 . Our estimate of the power-law component in the spectrum of 
SNR 1987A may be translated to a 2a upper limit for the luminosity of synchrotron emission of 
~ 3 x 10 35 erg s _1 in the 20-60 keV band. This value can be also considered as the limit for the 
luminosity of the stellar remnant (pulsar or black hole) harboured inside the ejecta. 
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Figure 1. Hard X-ray images indi- 
cating the detection of 4 Ti emission 
lines from SNR1987A. a-c, Maps of the 
signal-to-noise ratio (S/N) of the 1°5 x 
1°5 sky region around SNR1987A ac- 
cumulated in three energy bands with 
the IBIS/ISGRI telescope on board INTE- 
GRAL during observations in 2003-201 1 
(~ 6.0 Ms of real exposure or ~ 4.2 Ms 
of dead-time-corrected exposure): 48- 
65 keV (a); 65-82 keV (b); 82-99 keV 
(c); The maps were reconstructed us- 
ing standard technique^ with contours 
given at S/N levels of 2.7, 3.3, 3.9, 4.5, 
5.4 and 6.3. Two well-known sources, 
PSR B0540-69 and LMCX-1, are seen 
bright in all three images, but SNR 1987A 
is confidently detected only in b, in the 
band that contains the 67.9- and 78.4- 
keV direct-escape lines of radioactive 44 Ti 
decaying inside the ejecta. 
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Figure 2. Hard X-ray spectrum of SNR1987A measured with IBIS/ISGRI in 2003-2011. a, 

The observed (count) spectrum together with the simulated response curves; b, the unfolded (pho- 
ton) spectrum. In each panel, open circles show the spectra, dotted curves show the fit with two 
Gaussian lines at 67.9 and 78.4 keV of the radioactive decay of 44 Ti and solid curves show a similar 
fit which includes an additional power-law continuum with photon index a = 2.1. Each data point 
(open circles) represents the result of flux measurement at the position of SNR 1 987A in the image 
corresponding to the given energy band. It is expected that 44 Ti is produced inside the supernova 
core, which is expanding with velocity w 44 < 1, 700 km s _1 ; thus, the internal width of the lines (that 
is, their width before smoothing due to the finite energy resolution of the instrument) is unlikely 
to exceed AE ~ 0.4 keV. We kept the width of the lines used for fitting fixed at 2.8 keV FWHM, 
which is half of the energy resolution of ISGRI (such choice does not affect flux measurements in 
the internally narrow lines, thereby making it possible to work with the discrete response matrix of 
ISGRI). The upper limits are 1.7<r (90% confidence), error bars on the other points are la. 
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Figure 3. 7-ray spectrum of SNR 1 987A measured with SPI in 2003-201 1 . The spectrum was 
reconstructed using standard techniques 28 near the highest-energy (1 ,157 keV) direct-escape line 
of 44 Ti. The solid curve shows the Gaussian line with the flux at its 1.7a upper limit. Its centroid 
energy was redshifted by 1 .1 keV to take into account the heliocentric velocity of SNR 1 987A (285 
km s -1 ; ref. [29); its FWHM was taken to be 7.0 keV (internal width ~ 1,157 keV x v u /c ~ 
6.5 keV, corrected for the energy resolution of SPI, namely ~ 3 keV FWHM at 1.2 MeV; here 
u 44 < 1700 km s _1 is the expansion velocity of the ejecta layer containing 44 Ti, and c is the velocity 
of light). The similar Gaussian line with total flux corresponding to the 44 Ti mass measured by 
IBIS/ISGRI in the hard-X-ray lines is shown by a dotted curve. Error bars are la. 



